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Abstract 

The powers of matrices with Stirling number-coefficients are investigated. It is revealed 
that the elements of these matrices have a number of properties of the ordinary Stirling 
numbers. Moreover, "higher order" Bell, Fubini and Eulerian numbers can be defined. 
Hence we give a new interpretation for E. T. Bell's iterated exponential integers. In 
addition, it is worth to note that these numbers appear in combinatorial physics, in the 
problem of the normal ordering of quantum field theoretical operators. 
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Stirling cycle numbers (also known as Stirling number of the first kind [1]) count 
the number of permutations with given cycles. Concretely, |s„ „,| (also denoted as M or 
c(n, to)) gives the number of n-paritions with exactly to cycles (0 < m < n). The signed 
Stirling cycle number is {-l)"~ m \s nm \ =: s n m . The exponential generating function for 
the latter is 
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Stirling set number (also known as Stirling number of the second kind 11'!]), denoted 
by S„ tm (or {"J), enumerates the partitions of an n-set into to subsets (0 < to < «). The 
exponential generating function for them is 
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We enumerate the main identities involving these numbers. Let (x 
1) • • • (x — n + 1) (with (x)o = 1), then the next matrix identities are valid [ 1 
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So the Stirling matrices are inverses of each other. 

Applying again the matrix on (0, we get the Bell polynomials J2l: 



S o,o 
S i,o S 1,1 

S n ,0 $n,\ 









That is, 



S n,n-l S n n ) 



Bn(,X) — ^ ^ S n,m,X 



\ (x°} 




(Bo(x)\ 






Bi(x) 


v" 

/ \x ) 




,B„(x), 



(2) 



(3) 



(4) 



As one can see by the definition, B„ :- B„(l) gives all the possible partitions of a 
set with n element. These numbers are called Bell numbers |2J|. Among others, we 
point out that powers of the Stirling matrices present new numbers which satisfies all 
the relevant properties of the Bell numbers. It seems that these numbers were not 
investigated except by E. T. Bell himself [3]. 

In all what follows, we would like to discuss the Stirling matrices on a unified base 
as far as possible. So we introduce the notation 
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and the infinite lower triangular matrices 
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S n l m , with the notation (0). 



and S' 1 , which is the inverse of S. That is, (S^ 1 ),^, 

The formal p-th power of these matrices will be denoted as S p and the inverse is 
S~ p , while their («, m)-th entries are Sf l m and S n p m , respectively. 



1. Basic properties of the entries of S ±p 

First, to discuss the generating functions on the most simple way, we introduce the 
functions 

cr'O) := e x - 1, cr p (x) := cr^ 1 (cr 1 (*)), 
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and 



a-~\x) := ln(l + x), (r~ p {x) := o-^'V^W), 



here p > 1. Further, let <x°(x) = x. 
Then we have the general 

Proposition 1. The generating function of the m-th column ofS p (p e ZJ is 
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Proof. For p £ {— 1,0, 1} this is known or trivial. Let p > 1. Then 
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by induction. The case of p < — 1 is totally similar. 

2. Higher order Bell numbers 

Next we show that the well known identity JH 
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can be generalized and thus the "higher order" Bell numbers and polynomials defined 
below appear naturally. To see this, we cite Bell's definition [3] for the so-called "iter- 
ated exponential integers" - mutatis mutandis. 

Definition 2 (E. T. Bell - 1938). The higher order Bell polynomials ( or iterated expo- 
nential polynomials) can be defined via their exponential generating functions. Let 

E p {t) := 

where p > 0. Then we can define the polynomials 

B p n (x):=nl[f](E p (t)y, (7) 
and the "higher order" Bell numbers 

B p := B p (l). 

It can be seen that B®(x) = x". 

We give the following table of the numbers B p . In the paper of Ginsburg [4] there 
is a misprint. B^ was given incorrectly by him. 
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We remark that Bell used the notation and indices instead of B p n but the latter 
one fits in with our approach. He did not use the polynomials but only the numbers and 
gave several properties for B p . On the other hand, Ginsburg used the letter E also for 
the generating functions of the iterated exponentials but he applied lower indices. 

It should also be noted that these numbers appear in a paper concerning combina- 
torial physics: normal ordering of boson annihilation and creation operators Oh- 

With these considerations, we present the next 

Proposition 3. For p > we have the more general form of (0 as 
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and for p > 1 we get the generalization of 

si m = ±f J h(-ir- k B*-\k). 

Proof. We consider the series 
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we can continue as 
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Definition © and (fj) gives the wanted identity. 
To see the second one, 
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Comparing the coefficients, the result follows from PropositionQ] □ 
In his paper, Bell cited Dobinski's result [6] 
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and wrote that "there are similar (but more complicated) summations involving g„ , 
m > 1". Without cited references with respect to the mentioned formulas, we do not 
know what would be these complicated sums, but we point out that beside Dobinski's 
nice formula, the same and simple summation exist for higher order Bell numbers. 
What is more, Cesaro's integral representation formula [7] can also be carried to the 
more general setting. 

Proposition 4. For p > 1 we have 

1 A x k B p -\k) 
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Proof. The case p = 1 is the result of Dobihski. Let p > 1 . By the definition of 

B p n (x), 

y B p (x)- = (E p (t)Y = e** m = = 

_ 1 y jV^f _ 1 y X k [EP-\t)f 

~ e x ^-i k\ e x ^-i k\ 
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And the integral representation of higher order Bell numbers is as follows. 
Proposition 5. For p > we have 

2«' C" 

B p (x) = — -Im I exp(xo- p (e")) sin(nt)dt. 



Proof. The case p — comes from the orthogonality of sines on [0, n] and the 
Taylor expansion of exp and the DeMoivre's formula, p > is the formula of Cesaro, 
so let p > 1 . 
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3. Higher order Fubini numbers 



An other direction of generalizing old identities arises if we consider the "ordered 
version" of Bell numbers. 

The ordered Bell number (also known as Fubini number (H @, [l(J) F n , enumer- 
ates the ordered partitions of an n-set. Thus, we can define F„'s as the values of the 
polynomials 
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at x — 1 . The generating function is iTToL [Till 
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Nice identities are valid for these numbers and polynomials 
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where A„ ^ is the Eulerian number with parameters n and k and counts the permutations 
consisting k rises on n elements (see fljj] ) . 

In what follows, we show that these relations remain true in the more general con- 
text. 

Definition 6. Let p > 1, and let 
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We call these polynomials as higher order Fubini polynomials, and F^(\) —: F„ as 
higher order Fubini numbers. 

It is very interesting that these numbers also appear in d^. 
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Proposition 7. With the above definition, the Fubini numbers of higher order can be 
represented as 

B P n (k) 
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2 k + l ■ 



Proof. Rewrite Proposition^ as 

I! 

Z Sn,m x — Bn( x )- 
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If we apply Proposition!!] multiply both sides with l/e x and integrate on [0, +°o[ with 
respect to x, we get 
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The integral on the left hand side is ml, and 

r° _ /k! 

hence the result follows. □ 
It is not hard to deduce the exponential generating functions for these numbers. 

Proposition 8. We have for all p > 1 



00 n 1 

M=0 



JK7*(f) 

Proof. Case /? = 1 is known (see © and the references before the formula), so let 
p > 1 be fixed. We define a linear operator L on M[[x]] as 

L((x) ffl ) := mix"'. 

(This defines L entirely, because the sequence ((x) m )°°_ is a base in K[[x]].) Then, we 
apply L to ©: 



L(Br'(x)) = ^S^ZXMm) = ^S£ m m!x ffl = F p „(x). 

m=0 m=0 

This gives 
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Now we introduce the temporary variable v :- E p (t) - 1 (t can be considered as a 
fixed real number). Whence 

L(^'(0r) = K(l +v« = L Z^ - Yj^n\ = Z = 
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1 1 



1 - xv 1 - x(£ , ''- I (0 - 1) 1 - xo-P{t) ' 
This proves the proposition. 
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4. Higher order Eulerian numbers 



Now we show that how can one construct higher order Eulerian numbers. Let us 
recall ( fTTT i. and write it in the form 



x° 



:= M 



Here the (n, k)-th entry of J{ is A n ^, the Eulerian number with parameters n and k. 
According to Definition [6] and identity ©, it seems to be fruitful to give a similar 

Definition 9. Let ,9v p ' be the ( unique ) matrix (p > 1 ), for which 
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In addition, we introduce the polynomials 
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We shall point out that this definition is the proper one, since it gives back all the 
existing identities connecting Fubini, Stirling and Eulerian numbers and polynomials 
(PropositionfTOl). 

It is important to notice that 3^ is not equal to the p-th power of 3K. This is why 
we use the brackets. 

Because the identity (fTTT) was generalized on an other direction in a former paper 
of the author in [13], all the proofs of the next Proposition can be deduced verbatim. 
These identities are straight generalizations of the ordinary case (p = 1). 

Proposition 10. The next identities are valid 

k=0 N ' 

n 

F p n {x) = £V,(x+l)V-*, 

k=0 

F p n {x) = SA'l—). 
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Corollary 11. The exponential generating function for the polynomial A p n (x) is 

YK( X £ = * = *" 1 



re! \-A-o-P{{ x -Y) t ) x- EP- l ((x-l)t) 

5. Application to the eigensequence of the Stirling transform 

In this section we present a formula for the eigensequence of the Stirling transform 
involving S* m , First, let us introduce the necessary notions. Let a sequence a„ be 
given. Then the Stirling transform b„ of a„ is defined by 



b„ = ^a k S n , k (re = 0,1,...). 



k=0 



We consider the eigensequence of the Stirling transform 11411 . that is, C„ is defined as 

n 

C n +\ — j CkS n ,k- 
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(The sequence begins with 0, 1, 1, 2, 6, 26, 152, 1 144, ) With our notation, 
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It is known [14] that the exponential generating function of C„ satisfies the equation 

f'(x) = f(e x - 1) + 1, 
or, with the notations have introduced, 

f(x) = f(o-\x)) + l. 
Then it can be readily checked that 

/"(*) = e x [f{o-\x)) + l] 

We see that 

_2/VY>* °° °° 



(12) 



i CO 



k=0 k=0 n=0 

The coefficients of x"/n\ are the elements of the vector 
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(truncated after n). The final consideration is that 

/7=0 «=0 «=() 

Hence < fT2b implies the identity 

jfc=0 v 1 

or, rephrased applying the matrix formalism, 
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In addition, we prove the next 
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Proof. Let us introduce the operator £ as 
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Recalling (0, we get 



So, by [2, p. 114], 

C„ +2 = £B„ +l (x) = £(xB,,(x) + xB'Jx)) = £(xB„(x)) + £(xB' n (x)) 
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J]s„, k x k+[ +£j]s n , k kx k =Y j S n , k £(x k+l ) + Y j S n , k k£(x k ), 



\n=0 ' 

and the result follows 
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□ 



References 



[1] R. L. Graham, D. E. Knuth, O. Patashnik, Concrete Mathematics, 
Addison- Wesley, 1994. 

[2] M. Aigner, Combinatorial Theory, reprint of the 1979 edition, Springer- 
Verlag, Berlin Heidelberg, 1997 

[3] E. T. Bell, The iterated exponential integers, Ann. Math. 39(3) (1938), 
539-557. 

[4] J. Ginsburg, Iterated exponentials, Scripta Math. 1 1 (1945), 340-353. 

[5] K. A. Penson, P. Blasiak, G. Duchamp, A. Horzela, A. I. Solomon, Hi- 
erarchical Dobiriski-type relations via substitution and the moment prob- 
lem, J.Phys. A: Math.Gen. 37 3475-3487 (2004). 

[6] G. Dobihski, Summirung der Reihe 2« m /«! fur m = 1,2,3,4,5,..., 
Grunert Archiv (Arch, fur Mat. und Physik) 61 (1877), 333-336. 

[7] M. E. Cesaro, Sur une equation aux differences melees, Nouv. Ann. 
Math. 4 (1885), 36-40. 

[8] O. A. Gross, Preferential arrangements, Amer. Math. Monthly, 69 (1962), 
4-8. 

[9] R. D. James, The factors of a square-free integer, Canad. Math. Bull. 1 1 
(1968), 733-735. 

[10] S. M. Tanny, On some numbers related to the Bell numbers, Canad. Math. 
Bull. 17 (1975), 733-738. 

[11] R. Sprugnoli, Riordan arrays and combinatorial sums, Discrete Math. 
132(1994), 267-290. 



[12] H. S. Wilf, Generatingfunctionology, Academic Press, 1994. 



11 



[13] I. Mezo, On a question of Benoumhani related to Whitney numbers, 
Seminaire Lotharingien de Combinatoire (submitted). 

[14] M. Bernstein, N. J. A. Sloane, Some canonical sequences of integers, 
Linear Alg. Applications, 226-228 (1995), 57-72; erratum 320 (2000), 
210. 



12 



